Superconductor-insulator transition from phase fluctuations in layered 

superconducting films. 
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Phase fluctuations in finite thickness layered superconducting films are studied theoretically. The 
model is shown to reproduce much of the experimental findings observed in various systems, mainly 
a superconductor-insulator transition and an inverse dependence of T c on thickness, T c ~ d . The 
relevance to disordered thin films is discussed. 



The superconducting (SC) state and the insulating 
state are two extremes of the electronic conduction spec- 
trum. It is thus rather surprising that a direct tran- 
sition between them exists. However, various SC sys- 
tems such as amorphous and granular films of various 
compounds exhibit just that a SC to insulator transition 
(SIT) tuned by an external parameter, typically an exter- 
nal magnetic field of the nominal film thickness [1]. The 
details of the SIT have been a challenge to the community 
for more than two decades and still many open question 
remain, especially regarding the thickness-induced tran- 
sition. The standard picture is that as the film thickness 
is reduced the amount of disorder increases, eventually 
causing enough fluctuations that cooper-pairs cannot sur- 
vive, and thus the SC state is destroyed. The transition 
is viewed as a quantum phase transition [2, 3]. 

Here we introduce a very simple model which exhibits 
an apparent SIT and also captures other experimental 
finding, such as the dependence of T c on thickness. While 
not directly related to any specific experimental system, 
the phenomenological resemblance between the model 
and various experiments suggests that the mechanism we 
propose for the SIT should not be ruled out, if only from 
the principle of Occam's razor. 

Consider a system composed of stacks of truly 2D SC 
films. Being two dimensional, the action for the SC phase 
of each film (if they are independent) is 



(i) 



where <j> is the phase of the order parameter and K is the 
phase stiffness, which is proportional to the supcrfluid 
density, and SVortices describes the proliferation of vor- 
tices [4, 5]. This action is known to lead to a Kostcrlitz- 
Thoulcss-Bcrezinskii (KTB) transition at a temperature 
T c ~ . We assume that the only role of disorder is 
to determine the magnitude of the phase stiffness (the 
stronger the disorder the smaller the stiffness), and we 
do not assume any electronic inhomogeneity in the sys- 
tem (which may arise with an applied perpendicular field 
[6]). Since T c is determined by K we omit the Svorticcs 
term from the following actions and consider the effect 
of thickness on K . 



The next step is to assume some form of coupling be- 
tween the phases of the layers. The simplest form which, 
as we show below, can be solved exactly, is 
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^(V^(r)) 2 +J(0 l (r)-^_ 1 (r)) 2 + 



-JiV0i(r) ■ V<fc_i(r)] 
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where 4>i ( r ) is the phase of the order parameter in the i-th 
layer, J describes a a discrete generalization of the gradi- 
ent term to the z-direction, and J\ is a Bio-Savart type of 
interaction, which (as we show below) emerges naturally 
from the procedure we use. In an isotropic system one 
would expect that J /a 2 ~ Kq where a is the distance be- 
tween layers (in such a way that in a continuum limit this 
term restores the gradient term in the z-dircction), and 
that J\ be much smaller than Kq. Similar actions were 
previously studied using a renormalization group (RG) 
approach in the limit of infinite thickness [7, 8], where 
various crossovers were studied as a function of, e.g. the 
interlayer coupling J. 

The action of Eq. (2) is quadratic, and therefore it is 
possible to solve for an infinite system. However, we are 
interested in the effect of finite thickness, i.e. a finite 
number of layers, which can be obtained as follows. Con- 
sider the action of one of the layers, say with an even 
index 2i. Fourier transforming to momentum space, it is 
given by 
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where (f>^k is the k— th Fourier component of (j>2i{r), and 
\<Pk\ 2 = <Pk<P-k- It is now straight- forward to integrate 
out 4>2i,k and obtain an effective action (keeping only 



2 



second order in k) for the odd layers 
K 



(V(foi+l)~ + J{<t>2i+l - <t>2i-l)' 



-JiV02i+iV0 : 
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with K = K + (K /2 - 2 Ji), Ji = Ji - A'o/4, J = J/2. 
We now repeat the above procedure for the new cou- 
plings. This real-space RG procedure yields the following 
recursion relations, 



■) , (5) 



= ^_ 2J («) 



J(n+1) = j(n) /2 ^(n+1) = (J (n) 



with jW = J, A^ 1 ) = A , = Ji. The film thickness 
d is related to the RG index n via 2 n = - = 5. 

a 

The recursion equations above can be easily solved. 
For J we find simply = 2~™ = (<5)~ , i.e. it is 

decoupled from K and J±. The interesting consequence 
of the recursion relation is in K^ n \ which is given by 
K (n) = | 2 -™(A + 4Ji) + ±2 n (A - 2Ji). For J 1 < A 
(but also for the general case) , and in the limit of a many 
layers layer, we find 



K{d) ~ -A 5 



(0) 



Before we consider the consequences of Eq. (6), Let 
us first recap the phenomenology observed in typical ex- 
periments which exhibit the thickness-induced SIT [1, 9]. 
In the experiments, the resistance is measured as a func- 
tion of temperature for different film thickness. For thick 
films, thicker than some critical thickness d c , there is a SC 
transition (i.e. resistance becomes vanishingly small) oc- 
curring at some critical temperature T c which is thickness 
dependent. For films thinner than d c the resistance seems 
to increase with decreasing temperature, never reaching 
a zero-resistance state. Other than that, most details of 
this transition (i.e. the resistance at the transition, the 
existence of a temperature-independent seperatrix, criti- 
cal exponents etc.), are non-universal and depend on the 
particular system examined [9, 10]. Different Chemical 
composition, growth and annealing process affect these 
parameters.. 

Assume now that the experiments are performed down 
to a certain minimal temperature T m ; n (this is always 
correct). Since such an experiment does not probe what 
happens below T m i n , Eq. (6) means that for films thicker 
than 
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the temperatures of the experiment will drop below the 
film's T c , and it will be considered as a SC film, while 



thinner films would appear insulating (see Fig. 1(a)) . 
Above the critical temperature, the resistance of films 
above and below d c will be qualitatively similar (see the 
inset of Fig. 1 of Rcf. [11], where the resistances arc plot- 
ted on a linear scale) , and is determined by a complicated 
interplay of vortices, disorder and finite size effects [12]. 
A direct prediction of this argument is that for exper- 
iments with increasingly lower temperatures, the resis- 
tance of films just below d c will eventually start decreas- 
ing and convert to a SC state. Thus, the critical thickness 
will increase linearly with T m i n and samples which are 
currently thought to be insulating will turn out to be su- 
perconducting. This implies that the resistance will drop 
orders of magnitude within a small temperature range in 
the vicinity of the critical temperature [13]. 

The nature of the true T = ground state of these 
films and the possibility of a T = insulating state due 
to quantum fluctuations are still open questions. A pos- 
sible scenario is that quantum fluctuations give rise to 
an effective minimal temperature Tq to which the sys- 
tem can reach, which plays the role of T m j n as in the 
discussion above. This possibility will be explored in the 
future. 

We point that although we have presented a model 
which gives the result T c ~ d analytically, one can also 
understand this from the experimental data and the fact 
that in a disordered superconductor the superfluid stiff- 
ness is proportional to the mean-free path. The experi- 
mental data points that the thickness is simply control- 
ling the mean-free path, implying that the stiffness (and 
hence the critical temperature) is proportional to the 
thickness. Thus, the SIT mechanism wc propose may be 
relevant even in the general case which cannot be mapped 
to a layer model. 

From Eq. (6) one can also determine the thickness 
dependence of T c (see Fig. 2 of Rcf. [9]). Near the 
transition An, is approximately given by [12] Kq(T) ss 
An(0)(l — jr-^), where T C: n is the mean- field transition 
temperature. Inserting this into Eq. (6) gives an equation 
for T c for a given thickness d, A (0)(1 — = 6T c /ir . 

Solving this yields 

6T, 



for i > 1. Such a dependence have been observed long 
ago [9] and can be extracted from SIT data from various 
systems, e.g. Bi films [14], Ge and Xe films [15] and even 
high-Tc films [16] and seems quite universal. We are, 
however, unaware of any prior theoretical discussion for 
the origin of this dependence. 

From the above considerations it is simple to under- 
stand the presence of a the thermal dimensional crossover 
[7, 8]. The system is characterized by two transition 
temperatures, T c ^ (quasi-2D KTB transition) and the 
mean- field transition T c $, both decreasing with decreas- 
ing thickness (Fig. 1(b)). As long as T c ^ < T c .n the 
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FIG. 1: (a) The temperature-thickness phase diagram, 
demonstrating the consequence of a critical temperature 
which increases with thickness. In the experiments for each 
measurement the thickness is set and the temperature is low- 
ered down to a finite temperature T m i n . Consequently, for 
films thicker than d c (Eq. (7)) the temperature will cross T c 
and a SIT will take place. For films thinner than d c the 
temperature does not reach its critical value and the sam- 
ple remains in the insulating state, (b) Dependence of T c ^ 
and T c fi on thickness, indicating the crossover from 2D to 3D 
behavior. 



system will have a two-dimensional character (with re- 
sistance affected by vortices), but for thick enough films 
the transition will be a regular BCS mean-field type. The 
transport properties at the crossover region are of inter- 
est, and will be a subject of future study. 

To summarize, we have introduced a model for layered 
thin superconducting films. The main result of the model 
is the observation that the phase-stiffness, which is pro- 
portional to the critical temperature, is proportional to 
the film thickness. A direct consequence of this fact is 
a possible observation of an apparent SIT, and a depen- 
dence of T c on film thickness given by Eq. 8. 

We note that the model may be relevant to other SC 
systems, such as High-T c films (where recent progress in 
fabricating thin films was achieved [17]). These materials 
are naturally layered, and thus the theory presented here 
should adequately describe them. However, owing to the 
high bulk T c , they do not exhibit a T-SIT and remain SC 
down to very thin samples. We predict that a detailed 
thickness-dependence study in the extremely under- or 
over-doped regions should exhibit the same phenomenol- 



ogy as the T-SIT in disordered thin films. 
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